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Discrete Laplace-Beltrami Operator Determines 
Discrete Riemannian Metric 

Xianfeng David Gu^, Ren Guo*, Feng Luo*, Wei Zeng § 

Abstract 

The Laplace-Beltrami operator of a smooth Riemannian manifold is determined by the Riemannian metric. 
Conversely, the heat kernel constructed from its eigenvalues and eigenfunctions determines the Riemannian metric. 
This work proves the analogy on Euclidean polyhedral surfaces (triangle meshes), that the discrete Laplace-Beltrami 
operator and the discrete Riemannian metric (unique up to a scaling) are mutually determined by each other. 

Given an Euclidean polyhedral surface, its Riemannian metric is represented as edge lengths, satisfying triangle 
inequalities on all faces. The Laplace-Beltrami operator is formulated using the cotangent formula, where the edge 
weight is defined as the sum of the cotangent of angles against the edge. We prove that the edge lengths can be 
determined by the edge weights unique up to a scaling using the variational approach. 

First, we show that the space of all possible metrics of a polyhedral surface is convex. Then, we construct a 
special energy defined on the metric space, such that the gradient of the energy equals to the edge weights. Third, 
we show the Hessian matrix of the energy is positive definite, restricted on the tangent space of the metric space, 
therefore the energy is convex. Finally, by the fact that the parameter on a convex domain and the gradient of a 
convex function defined on the domain have one-to-one correspondence, we show the edge weights determines the 
polyhedral metric unique up to a scaling. 

The constructive proof leads to a computational algorithm that finds the unique metric on a topological triangle 
mesh from a discrete Laplace-Beltrami operator matrix. 

I. INTRODUCTION 

Laplace-Beltrami operator plays a fundamental role in Riemannian geometry [10]. Discrete Laplace- 
Beltrami operators on triangulated surface meshes span the entire spectrum of geometry processing 
applications, including mesh parameterization, segmentation, reconstruction, compression, re-meshing and 
so on [5], [9], [14]. Laplace-Beltrami operator is determined by the Riemannian metric. The heat kernel 
can be constructed from the eigenvalues and eigenfunctions of the Laplace-Beltrami operator, conversely, 
it fully determines the Riemannian metric (unique up to a scaling). In this work, we prove the discrete 
analogy to this fundamental fact, that the discrete Laplace-Beltrami operator and the discrete Riemannian 
metric are mutually determined by each other. 

Related Works In real applications, a smooth metric surface is usually represented as a triangulated 
mesh. The manifold heat kernel is estimated from the discrete Laplace operator. The most well-known and 
widely-used discrete formulation of Laplace operator over triangulated meshes is the so-called cotangent 
scheme, which was originally introduced in [3], [7]. Xu [13] proposed several simple discretization schemes 
of Laplace operators over triangulated surfaces, and established the theoretical analysis on convergence. 
Wardetzky et al. [12] proved the theoretical limitation that the discrete Laplacians cannot satisfy all natural 
properties, thus, explained the diversity of existing discrete Laplace operators. A family of operations were 
presented by extending more natural properties into the existing operators. Reuter et al. [8] computed 
a discrete Laplace operator using the finite element method, and exploited the isometry invariance of 
the Laplace operator as shape fingerprint for object comparison. Belkin et al. [1] proposed the first 
discrete Laplacian that pointwise converges to the true Laplacian as the input mesh approximates a smooth 
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manifold better. Tamal et al. [2] employed this mesh Laplacian and provided the first convergence to 
relate the discrete spectrum with the true spectrum, and studied the stability and robustness of the discrete 
approximation of Laplace spectra. The eigenfunctions of Laplace-Beltrami operator have been applied for 
global intrinsic symmetry detection in [6]. Heat Kernel Signature was proposed in [11], which is concise 
and characterizes the shape up to isometry. 

Our Results In this work, we prove that the discrete Laplace-Beltrami operator based on the cotangent 
scheme [3], [7] is determined by the discrete Riemannian metric, and also determines the metric unique up 
to a scaling. The proof is using the variational approach, which leads to a practical algorithm to compute 
a Riemannian metric from a prescribed Laplace-Beltrami operator. 

Paper Outline In Section II, we briefly overview the fundamental theorem of smooth heat kernel 
and our theoretical claims of discrete case. We clarify the simplest case, one triangle mesh in Section III 
first; then turn to the more general Euclidean polyhedral surfaces in Section IV. Finally, in Section V, we 
present a variational algorithm to compute the unique Riemannian metric from from a Laplace-Beltrami 
matrix. The numerical experiments on different topological triangle meshes support the theoretic results. 

II. Preliminaries and Proof Overview 

A. Smooth Case 

Suppose (Af,g) is a complete Riemannian manifold, g is the Riemannian metric. A is the Laplace- 
Beltrami operator. The eigenvalues {A„} and eigenfunctions {§ n } of A are 

A(j) n = -Xntym 

where <j>„ is normalized to be orthonormal in L 2 (M). The spectrum is given by 

= Xq < Ai < A2 < • • • , An — > °°. 
Then there is a heat kernel K(x,y,t) E C°°(M x M xR+), such that 

00 

K(x,y,t)= J £e- lnt Ux)Uy)- 

n=0 

Heat kernel reflects all the information of the Riemannian metric g. The details of the following theorem 
can be found in [11]. 

Theorem 2.1: Let / : (Mi,gi) — > (Af2,g2) be a diffeomorphism between two Riemannian manifolds. If 
/ is an isometry, then 

K 1 (x,y,t) = K 2 (f(x),f(y),t), Vx,yeM, t>0. (1) 
Conversely, if / is a surjective map, and Eqn. (1) holds, then / is an isometry. 

B. Discrete Case 

In this work, we focus on discrete surfaces, namely polyhedral surface. For example, a triangle mesh 
is piecewise linearly embedded in IR 3 . 

Definition 2.1 (Polyhedral Surface): An Euclidean polyhedral surface is a triple (S,T,d), where S is 
a closed surface, T is a triangulation of S and d is a metric on S whose restriction to each triangle is 
isometric to an Euclidean triangle. 

The well-known cotangent edge weight [3], [7] on an Euclidean polyhedral surface is defined as follows: 

Definition 2.2 (Cotangent Edge Weight): Suppose [v,-,v/] is a boundary edge of M, [vi,Vj] G dM, then 
[vi,Vj] is associated with one triangle [v,-, Vj, v*], the angle against [v,-,v/] at the vertex Vk is a, then the 
weight of [vi,Vj] is given by Wij = ^cota. Otherwise, if [vj,v/] is an interior edge, the two angles against 
it are a,/3, then the weight is w,j = ^(cota + cot/3). 

The discrete Laplace-Beltrami operator is constructed from the cotangent edge weight. 
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Definition 2.3 (Discrete Laplace Matrix): The discrete Laplace matrix L = (Ljj) for an Euclidean poly- 
hedral surface is given by 



Because L is symmetric, it can be decomposed as 

L = <Z>A<Z> T (2) 

where A = diag(Xo, Ai, • • • , A„), = Ao < Ai < A2 < • • • < A M , are the eigenvalues of L, and <!> = (0o|0i|02| • • • \<j> n ), 
L§i — A ? -0j, are the orthonormal eigenvectors, such that §J §j = 

Definition 2.4 (Discrete Heat Kernel): The discrete heat kernel is defined as follows: 

K(t) = <f>exp(-At)<f> T . (3) 
The Main Theorem, called Global Rigidity Theorem, in this work is as follows: 

Theorem 2.2: Suppose two Euclidean polyhedral surfaces (5,r,di) and (S,T,d2) are given, 

Li = L2, 

if and only if di and d2 differ by a scaling. 

Corollary 2.3: Suppose two Euclidean polyhedral surfaces (5,T,di) and (S,T, d.2) are given, 

K l (t)=K 2 (t),Vt>0, 

if and only if di and d2 differ by a scaling. 
Proof: Note that, 

dK(t) . 

Therefore, the discrete Laplace matrix and the discrete heat kernel mutually determine each other. ■ 

C. Proof Overview for Main Theorem 2.2 

The main idea for the proof is as follows. We fix the connectivity of the polyhedral surface (S,T). 
Suppose the edge set of (S,T) is sorted as E = {ei,e2,--- ,e m }, where m=\E\ number of edges, the face 
set is denoted as F. A triangle [v^v^vj G F is also denoted as {ij,k} G F. 

By definition, an Euclidean polyhedral metric on (5, T) is given by its edge length function d : E — >■ M + . 
We denote a metric as d = (d\,d2,- • ■ ,d m ), where di = d(ei) is the length of edge e\. Let 

E d (2) = {(di,d 2 ,d 3 )\di + dj > d k } 

be the space of all Euclidean triangles parameterized by the edge lengths, where {1,7,^} is a cyclic 
permutation of { 1 , 2, 3}. In this work, for convenience, we use u = U2, • • ■ , u m ) to represent the metric, 
where u k = \d\. 

Definition 2.5 (Admissible Metric Space): Given a triangulated surface (S,K), the admissible metric 
space is defined as 

m 

Q. u = {(mi,w 2 ,W3--- ,u m )\ Y,Uk = m,(V^i,y/uj,Vuk) G E d (2), V{/,7, k} e F}. 

k=\ 

We show that Q. u is a convex domain in W 71 . 

Definition 2.6 (Energy): An energy E : Q. u — > R is defined as: 

r(ui,u 2 -,u m ) m 

E(ui, U2--- , u m ) = / Y w k (n)dn k , (4) 

k=i 

where w^(ju) is the cotangent weight on the edge et determined by the metric /1. 



4 




Next we show this energy is convex in Lemma 3.5. According to the following lemma, the gradient of 
the energy V£(d) : Q ->■ R m 

VE : (ui,u 2 --- ,u m ) (wi,w 2 ,---w m ) 

is an embedding. Namely the metric is determined by the edge weight unique up to a scaling. 

Lemma 2.4: Suppose Q. C W is an open convex domain in W, E : D. — y M is a strictly convex function 
with positive definite Hessian matrix, then VE : Q. — y MP is a smooth embedding. 

Proof: If p ^ q in Q, let = (1 -t)p + tq e CI for all t G [0, 1]. Then f(t) = E(y(t)) : [0, 1] ->■ M 
is a strictly convex function, so that 

df{t) 



dt 



= V£| yW -(q-p). 



Because 



therefore 



^P = (q-p) r //| yW (q-p)>0, 



V£(p).(q-p)^V£(q).(q-p). 



This means V2s(p) 7^ V£(q), therefore VE 1 is injective. 

On the other hand, the Jacobi matrix of VE is the Hessian matrix of E, which is positive definite. It 
follows that VE : £1 — y MP is a smooth embedding. ■ 

From the discrete Laplace-Beltrami operator (Eqn. (2)) or the heat kernel (Eqn. (3)), we can compute 
all the cotangent edge weights, then because the edge weight determines the metric, we attain the Main 
Theorem 2.2. 

III. Euclidean Triangle 

In this section, we show the proof for the simplest case, a Euclidean triangle; in the next section, we 
generalize the proof to all types of triangle meshes. 

Given a triangle {i,j,k}, three corner angles denoted by 0,-, k }, three edge lengths denoted by 
{di.dj^d^}, as shown in Fig. 1. In this case, the problem is trivial. Given w^) = (cot0,-,cot0j,cot0 fc ), 

we can compute 0/, Ok) by taking the arctan function. Then the normalized edge lengths are given by 

3 

(di,dj,dk) = . _ ■ . _ ■ . - (sm0i,sm0j,sm0 k ). 
sin 0i + sin 0j + sin Ok 

Although this approach is direct and simple, it can not be generalized to more complicated polyhedral 
surfaces. In the following, we use a different approach, which can be generalized to all polyhedral surfaces. 
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(8) 



Lemma 3.1: Suppose an Euclidean triangle is with angles {6j,6j,6k} and edge lengths {di,dj,dk}, 
angles are treated as the functions of the edge lengths 6i(di,dj,dk), then 

d_0 L= d L 

ddi 2A K ' 

and 

-dd=-YA C0 ^ (6) 

where A is the area of the triangle. 

Proof: According to Euclidean cosine law, 

dl + d}-df 

we take derivative on both sides with respective to di, 

sinfldft = ~ 2d i 
1 ddi 2djdk 

ddi _ di _ di 
ddi djdk sin 0, 2A 

where A = jdjdk sin 0, is the area of the triangle. Similarly, 

d 9 9 9 d 

— (dy +d h -d i ) = —(2djd k cos 0,-) 

<90,- 

2J; = 2d k cos Qi — 2djdk sin 0*^-7- 
2A^— = <4cos 0, — dj = —di cos % 

We get 

<90, di cos 6k 
~ddj ~ 2A 

Lemma 3.2: In an Euclidean triangle, let m, = ^d 2 and w 7 - = then 

<9 cot 0, (9 cot 0j 

duj dui 

Proof: 

d cot di 1 <9cot0 ; 1 1 ddi 11 di cos 6k df cos 9k 



(9) 



<9m; <9Jj sin 2 Qi ddj dj sin 2 0,- 2A sin 2 0,- 2Adidj 

4R 2 cos 0* 



2A J/Jj 

where R is the radius of the circum circle of the triangle. The righthand side of Eqn. (10) is symmetric 
with respect to the indices i and j. ■ 

Corollary 3.3: The differential form 

CO = cot Oidui + cot Ojduj + cot Okduk (11) 
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Fig. 2. The geometric interpretation of the Hessian matrix. The in circle of the triangle is centered at O, with radius r. The perpendiculars 
rti, tij and are from the incenter of the triangle and orthogonal to the edge e/, ej and respectively. 



is a closed 1-form. 

Proof: By the above Lemma 3.2 regarding symmetry, 

d cot Oj dcotdi. , , .dcoidk dcotOj 
a(0 = ( — - — - )dui Aduj + ( — - - — -)duj Adu^ 

OUj OUj OUj OUfr 

d cot Oi dcoidk. . 

+ ( — =: ^ jdu^Adui 

dUk out 

= 0. 

■ 

Definition 3.1 (Admissible Metric Space): Let ui = jdf, the admissible metric space is defined as 

Q. u := {(ui,Uj,u k )\(<Jui,y/uj,y/uk) G E d (2), Ui + Uj + u k = 3} 
Lemma 3.4: The admissible metric space Q. u is a convex domain in M 3 . 

Proof: Suppose (ui,Uj,Uk) G Q u and G Q u , then from y/ui+ y/uj > y/u^, we get w,- + wy + 

l^fuflFj > Ufa. Define (uf,uj,u^) = X(ui,Uj,Uk) + (1 — A)(w ( -, where < A < 1. Then 

u t Uj = (Xut + (1 — X)Ui)(Xuj + (1 —X)Uj) 

= X. 2 UiUj + (1 — X) 2 Uitij + X(\ — X)(uiUj + UjUi) 
> X 2 UiUj + (1 — X) 2 Uiiij > + 2X{\ —X)yJ UiUjUiiij 
= (X y/UiUj + ( 1 - X)y/uiUj) 2 

It follows 



uf + 2\J ufu^ > X (ui + Uj + 2^/UiUj) + (1 — X ) ( w,- + Uj + 2 yj Uiuf) 
> Xuk+ (1 — X)Uk = u£ 

This shows (uf,uj,ul) G Q M . ■ 

Similarly, we define the edge weight space as follows. 

Definition 3.2 (Edge Weight Space): The edge weights of an Euclidean triangle form the edge weight 
space 

Cl e = {(cot0/,cot0y,cotQfc)|O< 0/,0/,fy < M; + 0/ + k = %}. 

Note that, 

1 - cot 0, cot 0/ 



cote^ = -cot(0 ; + J ) 



cot 0,- + cot 0/ 
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Lemma 3.5: The energy E 

Aui,Uj,u k ) 

E(ui,Ui,U)r) = / cotdjdTi + cotOjdTj + co\.O k dz k 

7(1,1,1) 



(12) 



is well defined on the admissible metric space £2 M and is convex. 

Proof: According to Corollary 3.3, the differential form is closed. Furthermore, the admissible metric 
space Q. u is a simply connected domain. The differential form is exact, therefore, the integration is path 
independent, and the energy function is well defined. 
Then we compute the Hessian matrix of the energy, 



H = 



2R 2 



1 


cos 8 k 


cos Oj 




djdj 


did k 


COS 8 k 


1 


cos 6i 


djdi 


d ) 


djd k 


cos Oj 


cos 6j 


1 


d k di 


d k dj 





2R 2 



(ViM) inu^j) (n»i\k) 
(ij,id (ij^k) 

(lk,li) ilk,!]) (lk,lk) 



As shown in Fig. 2, J,n ; - + djiij + d k n k = 0, 



n 



rdi ' ^ rJy ' 1,1 rJfc 



,1k 



where r is the radius of the incircle of the triangle. Suppose (xi,Xj,x k ) G K 3 is a vector in R 3 , then 



[ x i ; x j ) ^fc] 



inuni) inunj) inunk) 
dj,!]) (ij,ik) 

(lk,rii) (lk,lj) (ihlk) 









X J 




. Xk . 



\xiili+xjrij + x k ri k \\ 2 >0. 



If the result is zero, then (xi,Xj,x k ) = X(ui,Uj,u k ),X G IL That is the null space of the Hessian matrix. In 
the admissible metric space Q. u , Ui + Uj - + u k = C(C = 3), then dui + duj + du k — 0. If {dui,duj,du k ) belongs 
to the null space, then (diii,duj,du k ) = X(ui,Uj,u k ), therefore, A,(ui + uj + u k ) = 0. Because ui,Uj,u k are 
positive, A = 0. In summary, the energy on Q. u is convex. ■ 

Theorem 3.6: The mapping VE : Q. u — > Q.Q,(ui,Uj,u k ) — > (cot0 ; ,cot0j,cot%) is a diffeomorphism. 

Proof: The energy E(ui,Uj,u k ) is a convex function defined on the convex domain Q. u , according 
to Lemma 2.4, VE : (ui,Uj,u k ) — > (cot0,-,cot0j,cot0 fc ) is a diffeomorphism. ■ 

IV. Euclidean Polyhedral Surface 
In this section, we consider the whole polyhedral surface. 

A. Closed Surfaces 

Given a polyhedral surface (5, T,d), the admissible metric space and the edge weight have been defined 
in Section II-B respectively. 

Lemma 4.1: The admissible metric space Q. u is convex. 
Proof: For a triangle {ij,k} G F, define 

&u k ■= {(m,u j ,u k )\(^u i ,,/u],^u k ~) G£ d (2)}. 

Similar to the proof of Lemma 3.4, Cl l J k is convex. The admissible metric space for the mesh is 

m 

&u = H &u k f){( u hU2r-- ,Um)\ Y< U k = m }i 
{i,j,k}eF k=l 

the intersection Q. u is still convex. ■ 
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Definition 4.1 (Differential Form): The differential form ft) defined on £l u is the summation of the 
differential form on each face, 

m 

ft) = Mijk = lw i du U 

{i,j,k}eF (=1 

where (0^ is given in Eqn. (11) in Corollary 3.3. w ( - is the edge weight on e,-. 
Lemma 4.2: The differential form ft) is a closed 1-form. 
Proof: According to Corollary 3.3, 

dco= Y d(Dijk = 0. 

{i,j,k}eF 

U 

Lemma 4.3: The energy function 

r(ui,u 2 , -,«„) n 

E(u l ,u 2 ,--- ,u„) = L E ijk (ui,u 2 ,--- ,u n ) = l^Widui 
{i,j,k}eF AM,-,!) (=i 

is well defined and convex on Q. u , where is the energy on the face, defined in Eqn. (12). 

Proof: For each face {i,j,k} G F, the Hessian matrices of is semi-positive definite, therefore, 
the Hessian matrix of the total energy E is semi-positive definite. 

Similar to the proof of Lemma 3.5, the null space of the Hessian matrix H is 

kerH = {X(di,d 2 , • • • ,d n ),X e R}. 

The tangent space of Q. u at u = (wi,W2,- -- )"«) is denoted by TQ. u (u). Assume (dui, du 2 , ■ ■ ■ ,du n ) G 
T£l u {u), then from YJJLi = m > we g et L^=i du m = 0. Therefore, 

TQ. u (u)f]KerH = {0}, 

hence // is positive definite restricted on T£l u {u). So the total energy E is convex on Cl u . ■ 

Theorem 4.4: The mapping on a closed Euclidean polyhedral surface VE : Q. u — y R m , (ui,u 2 ,--- ,u n ) — > 
(wi,w 2l --- ,w n ) is a smooth embedding. 

Proof: The admissible metric space Cl u is convex as shown in Lemma 4.1, the total energy is convex 
as shown in Lemma 4.3. According to Lemma 2.4, VE is a smooth embedding. ■ 

B. Open Surfaces 

By the double covering technique [4], we can convert a polyhedral surface with boundaries to a closed 
surface. First, let (S,T) be a copy of (S,T), then we reverse the orientation of each face in M, and glue 
two surfaces S and S along their corresponding boundary edges, the resulting triangulated surface is a 
closed one. We get the following corollary 

Corollary 4.5: The mapping on an Euclidean polyhedral surface with boundaries VE : Q. u — y R m , (u\ , u 2 , 
■■■ ,u n )— y{w\,w 2 ,--- ,w n ) is a smooth embedding. 

Surely, the cotangent edge weights can be uniquely obtained from the discrete heat kernel. By combining 
Theorem 4.4 and Corollary 4.5, we obtain the major Theorem 2.2, Global Rigidity Theorem, of this work. 
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Genus 



Genus 1 



Genus 2 



Fig. 3. Euclidean polyhedral surfaces used in the experiments. 



V. Numerical Experiments 



From above theoretic deduction, we can design the algorithm to compute discrete metric with user 
prescribed edge weights. 

Problem 5.1: Let (S,T) be a triangulated surface, w(VPi,m>2, • •• ,VP n ) are the user prescribed edge 
weights. The problem is to find a discrete metric u = («i,W2,-"' ,u n ), such that this metric u induces 
the desired edge weight w. 

The algorithm is based on the following theorem. 

Theorem 5.2: Suppose (5, T) is a triangulated surface. If there exists an u G Q. u , which induces w, then 
u is the unique global minimum of the energy 



Proof: The gradient of the energy VE(u) = w — w, and since VE(u) = 0, therefore u is a critical 
point. The Hessian matrix of E(u) is positive definite, the domain Q. u is convex, therefore u is the unique 



In our numerical experiments, as shown in Fig. 3, we tested surfaces with different topologies, with 
different genus, with or without boundaries. All discrete polyhedral surfaces are triangle meshes scanned 
from real objects. Because the meshes are embedded in M 3 , they have induced Euclidean metric, which 
are used as the desired metric u. From the induced Euclidean metric, the desired edge weight w can 
be directly computed. Then we set the initial discrete metric to be the constant metric - ,1). By 

optimizing the energy in Eqn. (13), we can reach the global minimum, and recovered the desired metric, 
which differs from the induced Euclidean metric by a scaling. 



We conjecture that the Main Theorem 2.2 holds for arbitrary dimensional Euclidean polyhedral man- 
ifolds, that means discrete Laplace-Beltrami operator (or equivalently the discrete heat kernel) and the 
discrete metric for any dimensional Euclidean polyhedral manifold are mutually determined by each other. 
On the other hand, we will explore the possibility to establish the same theorem for different types of 
discrete Laplace-Beltrami operators. 




(13) 



global minimum of the energy. 



VI. Future Work 
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